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We modified the way in which the Universal Map is obtained in the regular dynamics to derive 
the Universal a-Family of Maps depending on a single parameter a > which is the order of the 
fractional derivative in the nonlinear fractional differential equation describing a system experiencing 
periodic kicks. We consider two particular a-families corresponding to the Standard and Logistic 
Maps. For fractional a < 2 in the area of parameter values of the transition through the period 
doubling cascade of bifurcations from regular to chaotic motion in regular dynamics corresponding 
fractional systems demonstrate a new type of attractors - cascade of bifurcations type trajectories. 

PACS numbers: 05.45.-a, 05.45.Pq, 87.23. Cc 



Fractional differential equations (FDE) are frequently 
used in science and engineering to describe systems with 
memory (see [JJ [2] and references there). Circuit ele- 
ments with memory: memristors, memcapacitors, and 
meminductors (see Wikipedia [5]) can be used to model 
such systems. Nonlinearity plays an important role in 
such systems [TJ. We'll call systems which can be de- 
scribed by the FDE fractional dynamical systems (FDS) 
even though the group property is not satisfied in this 
case and the values of the system variables depend on the 
whole history of the system's evolution. Because FDE are 
integro-differential equations and there are no high order 
numerical algorithms to simulate such equations, deriva- 
tion of the fractional maps (FM) is important for the 
investigation of the general properties of the FDS. First 
FM were derived from the FDE in [5-8J. First results 
of the investigation of the FM (see [S [7] H [TO] ) revealed 
new properties of the FDS. Cascade of bifurcations type 
trajectories (CBTT) are the most unusual features of the 
investigated FM. In CBTT cascade of bifurcations is not 
a result of the change in a system parameter but appears 
as the attracting single trajectory and is a new type of 
attractors. All previous investigations of the FM were 
done on the various forms of the fractional 2-dimensional 
Standard Map corresponding to the order 1 < a < 2 
of the fractional derivative. CBTT appeared in all in- 
vestigated FM. The role of the cascades of bifurcations 
in the transition from order to chaos in the regular dy- 
namics and their connection to the scaling properties of 
the corresponding systems are well investigated (see [IT]). 
Consideration of the origin and the necessary and suffi- 
cient conditions of the CBTT's existence requires further 
investigation of the FM, which includes development of 
the simple, if possible one-dimensional, FM. The best in- 
vestigated one-dimensional regular map is the ubiquitous 
Logistic Map. This map, and the maps with a < 1 in 
general, can't be derived in a way previously used ([TJIE]) 
to derive FM for a > 1 (for the discussion on the way in 
which the Logistic Map can be derived from a differential 



equation see [H]). 

To derive the equations of the Universal a-Family of 
Maps (UaFM) let's consider the following equation: 

t Q 00 , 

~ + G K (x(t-AT)) S(--(n + e))=0, (1) 

n— — 00 

where e > A > 0, a £ K, a > 0, in the limit e 0. 
The initial conditions should correspond to the type of 
the fractional derivative we are going to use. In the case 
a = 2, A = 0, and Gr{x) = KG[x) Eq. ([TJ corresponds 
to the equation which integration produces the regular 
Universal Map. Case A = 0, Gk{x) = KG(x), a > 1 
has been used to derive the fractional Universal Map. 
A 7^ is essential for the case a < 1 when x{t) is a 
function discontinued at the time of the kicks [5] [J5] and 
the use of the K as a parameter rather than a factor is 
necessary to extend the class of the considered maps to 
include the Logistic Map. Without losing the generality 
we assume T = 1. Case T 7^ 1 is considered in [12] and 
can be reduced to this case by rescaling the time variable. 

In the case of the Riemann-Liouville fractional deriva- 
tive Eq. ([TJ can be written as 

00 

D?x(t) + G K (x(t-A)) ]T <f(i-(n + e)) =0, (2) 
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where £ > A > 0, e ^ 0, < JV - 1< a < iV, a € I, 
N € Z, and the initial conditions (o-D" -fe a;)(0+) = Cfc, 
k = 1,...,N. The left-sided Riemann-Liouville fractional 
derivative D"x(t) defined for t > [TML3] as 

D?x(t) = D? ir a x{t) 

1 d n f* x(r)dT 
T(n-a)dt n J Q (t-T) a - n + 1 ' 1 ' 

where n — 1 < a < n, neZ, £)" = d n /dt n , and qI" is 
a fractional integral. 

This problem can be reduced [TJ [T5J [TTj] to the Volterra 
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integral equation of the second kind for t > 

N 

x{t) = E 
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which integration gives (t > 0) 

*w = E 
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where 0(i) is the Heaviside step function. In Eq. ([5| 
we took into account that boundness of x(t) at t = 
requires c^r = and x(0) = 0. 

With the introduction [5] = a Df~ N+1 x(t), 

pO)(t) = P t s p(^), s = 0,1, ...,7V- 2 Eq. ([5]) leads to 
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where s = 0,l,...,iV — 2. With the definitions a;„ = x(n) 
and y>n = (n) Eqs. ([5]) and (|6| in the limit e — > give 
for t=n+l the Riemann-Liouville UaFM (UaRLFM) 
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Map equations for the momentum defined in a usual way 

p{t) = D]x{t), P s {t) = D s tP (t), a = 0,1,..., JV-2, (9) 

and the discussion on the different ways of the defining 
momentum in the case of the Riemann-Liouville maps 
can be found in [12]. UaRLFM Eqs. Q ^ can be writ- 



ten in the much simpler form 
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where s = 0, 1, ...JV-2 and U Q fe (m) = m a - k -(m-l) a - k . 
For the integer a — N the UaFM converges to 
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with the Jacobian (JV x JV, N > 2) 
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where n and fc are the column and row numbers. The 
first column can be written as the sum of the column 
with one in the first row and the remaining zeros and the 
column which is equal to Gk(x) times the last column. 
Determinants of the corresponding matrices are 1 and 0, 
this is why the Jacobian is equal to one and the map is 
the N-dimensional volume preserving map. 

The integer UaFM's fixed points are p s Q = (s = 
0, ...,JV — 2) and Xo from G(xq) = 0. Their stability 
for TV > 1 is defined by the eigenvalues A of the Ja- 
cobian matrix J(x ,Pq, ...,p^~ 2 ). Polynomial P(A) = 
det[J(xo,Po, ■■■iPo ~ 2 ) ~ nas values P(0) = Ai x ... x 
Atv = 1 and P(l) = (— l) N Gk(xq), which means that for 
odd values of JV > 1 stability is possible only if Gk(xq) = 
0. For T = 2 points p s n+1 = -p s n (s = 0, N - 2) and 
G(x n +i) = —G(x n ). In the case N — 3 the only T = 2 
points are the periodic points. 

For Eq. ([I]) with the left-sided Caputo derivative [T5] 



D^x(r)d 



= rv V / n T J-n+i (n-Ka<n) (14) 

r(n - a) 7 (t - r) a "+ 1 

the initial conditions may be taken as (D k x)(0+) = 
bk, k = 0, ...,iV — 1. This problem is equivalent to the 
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Volterra integral equation of the second kind (t > 0) 
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With the introduction x( s )(i) = D^x(t) the Caputo 
UaFM (UaCFM) can be derived in the form [T] 
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where s = 0, 1, N — 1. 
Fractional maps Eqs. (10), (111, and (16) are maps 



with memory in which the next value of the map variables 
depends on all previous values. An increase in a leads to 
the increase in the dimension of the map and to the power 
law increase of the weights of the old states (the increased 
role of memory). Integer values of a correspond to the 
degenerate cases in which map equations can be written 
as maps with full memory |17j which are equivalent to the 
one step memory maps (for a discussion on the fractional 
maps as maps with memory see |12|h 

In the a = 2 case Eqs. ( 12 ) and ( 13 ) produce the Stan- 
dard Map if Gk{x) = Ksixx(x) and in the a = 1 case 
the Logistic Map results from Gk(x) = x — Kx(l — x). 
We'll call the UaFM Eqs. flO} and dTTh with G K (x) = 



Ksin(x) the Standard a-RL-Family of Maps (SaRLFM) 
and with Gk{x) — x — Kx(l — x) the Logistic a-RL- 
Family of Maps (LaRLFM); we'll call UaFM Eq. ^ 
with Gk{x) = K sm(x) the Standard a-Caputo-Family 
of Maps (SaCFM) and with G K (x) =x- Kx(l ~ x) the 
Logistic a-Caputo-Family of Maps (LaCFM). 

For a = the solution of Eq. ([I]) is identical zero. 
For a < 1 the UaRLFM Eq. ^ also produces identical 
zero for maps which satisfy G(0) = 0, which is true for 
the SaRLFM and LaRLFM. The a = 1 SaRLFM is a 
particular form of the Circle Map with zero driving phase 



x n+ \ — x n — K sm(x n ), (mod 2ir) 



(17) 



(see its bifurcation diagram Fig. Qa)). In this map the 
stable T = 4 sink appears at K w 3.445 and the tran- 
sition to chaos through the period doubling cascade of 
bifurcations occurs at K « 3.532 [T^j. In the a = 2 
Standard Map 



Pn+l = Pn - K sin X, x n+1 = x n + p n+1 



(18) 



T = 4 stable elliptic points appear at K w 6.59 and the 
period doubling cascade of bifurcations leads to the dis- 
appearance of the islands of stability in the chaotic sea at 
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FIG. 1. (a) The bifurcation diagram for ID Standard Map 



Eq. 17 (b) On K-a graph in SaRLFM and SaCFM CBTT 



exist in the band of the map's parameters ending at the cusp 
in the top right corner. Star marks the point at which the 
Standard Map's (a = 2) T = 2 elliptic points with x n +i = 
x n — tt and Pn+i = — Pn become unstable and bifurcate. 
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FIG. 2. A single CBTT in the SaRLFM. (a) One of the 
two branches of the CBTT for a = 1.1, K = 3.5. (b) A 
zoom of a small feature in an intermittent trajectory for a = 
1.95, K — 6.2. (c) An intermittent trajectory in phase space 
for a = 1.65, K — 4.5. (d) x of n for the case (c). 



K « 6.6344 PI|. In SaRLFM and SaCFM for 1 < a < 2 
CBTT exist in the band between two curves connecting 
the above-mentioned points Fig. [ljb) . The lower curve 
is calculated semi-analytically as the condition of the ap- 
pearance of T = 4 sink in the limit n — > oo. Both curves 
are confirmed by the large number of computer simula- 
tions [HI HO]- Within the CBTT band trajectories evolve 
from being very stable features which exist for the longest 
time we were running our codes, 500000 iterations, when 
a is close to one Fig. [2]ja) to being barely distinguishable 
and short-lived features when a is close to two Fig. [2jL) . 
For the intermediate values of a CBTT behave similar 
to the sticky trajectories in Hamiltonian dynamics: oc- 
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FIG. 3. Bifurcations and CBTT in SaCFM and LaCFM with 
< a < 1. (a)-(d) bifurcation diagrams obtained after per- 
forming 10 4 iterations on a single trajectory with xq — 0.1 for 
various values of K. (a) SaCFM with a = 0.5. (b) LaCFM 
with a = 0.5. (c) SaCFM with a = 0.05. (d) LaCFM with 
a = 0.1. (e) A CBTT in SaCFM with a = 0.01 and K = 276. 
(f) A CBTT in LaCFM with a = 0.1 and K = 22.7. 



casionally trajectories enter CBTT and then leave them 
entering the chaotic sea (Figs. [2|c),(d)). More on the 
properties of the a — 1 Standard Map and the consis- 
tency of the properties of the SaRLFM and SaCFM for 
1 < a < 2 can be found in [12]. 

There are no stable fixed points in the a = 3 Standard 
Map. For K 2 — 16 < 4p l2 < K 2 there exist two lines of 
the stable on the torus T = 2 points. For more on the 
preliminary results of the investigation of the SaRLFM 
and SaCFM for 2 < a < 3 see [12]. A different form of 
the 3D Standard Map has been recently introduced and 
investigated in 19 . SaRLFM and SaCFM for a > 2 are 
poorly investigated and 3D volume preserving maps are 
not fully investigated. 

With the corresponding Gr (x) UaCFM for < a < 1 
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FIG. 4. An inverse CBTT in the LaRLFM with a = 1.15, 
K — 3.45. 60000 iterations on a trajectory with xo = 0.01 
and po — 0.1. (a) Phase space, (b) x — n graph. 



produces SaCFM and LaCFM, which are one dimen- 
sional maps with the power law decreasing memory. The 
bifurcation diagrams for these maps are similar to the 
corresponding diagrams for the a = 1 case but are 
stretched along the parameter if-axis and the stretch- 
iness increases with the decrease in a Figs. [3]^a)-(d). In 
the area of the parameter values for which on the bifurca- 
tion diagram stable periodic T > 2 points exist individual 
trajectories are CBTT Figs. |);e),(f). 

The area preserving quadratic a = 2 Logistic Map 



Pn 



■ p n + Kx n (l-X n )- 



Xn+1 = X n +p n+1 (20) 



fc=0 



is similar to the map studied by Henon [20] (for a recent 
review see [3T] ). The map Eq. |20| ) has two fixed points: 
(0,0) stable for K € (-3,1) and ((K - 1)/K,0) stable 
for K € (1)5). At K = 5 as a result of a bifurcation 
a couple of T = 2 stable points appear. For more on 
the a = 2 and a = 3 quadratic Logistic Maps see [12]. 
3D quadratic volume preserving maps were investigated 
in [55]. Logistic families of maps for 1 < a < 3 are 
poorly investigated. Results of preliminary simulations 
show that CBTT are present in those maps. Fig.|4]shows 
an inverse CBTT in LaRLFM with a = 1.15, K = 3.45. 

Summary. The Universal a- Family of Maps introduced 
in this paper is the extension of the fractional Univer- 
sal Map, which allows consideration of the Logistic Map 
as its particular form. Using the Standard and Logistic 
Families of Maps we showed that the existence of the cas- 
cade of bifurcations type trajectories is a general property 
of the fractional dynamical systems. They appear for the 
parameter values corresponding to the transition through 
the period doubling cascade of bifurcations from regular 
to chaotic motion in the regular dynamics. Fig. [2] and 
Fig. [3] support our statement that with the increase in a, 
which represents the increase in the systems' dimension 
and memory (increase in the weights of the earlier states) , 
systems demonstrate more complex and chaotic behav- 
ior. Biological systems are systems with memory and 
the Fractional Logistic Map can serve as a basic model 
in population biology with memory. New types of mate- 
rials with memory, such as memristors, memcapacitors, 
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and meminductors, can be used to model fractional sys- 
tems to demonstrate the existence of the CBTT. a > 2 
Standard and Logistic Maps (including their integer vol- 
ume preserving forms) are topics of the ongoing research 
and their further investigation is necessary to demon- 
strate the consistency of the changes in the properties of 
the fractional systems with the change in a. 

The author expresses his gratitude to V. E. Tarasov for 
the useful remarks and to E. Hameiri and H. Weitzner 
for the opportunity to complete this work at the Courant 
Institute. 
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